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We calculate discrete beta functions corresponding to the two-lattice matching for the 2D 0{N) 
models and Dyson's hierarchical model. We describe and explain finite-size effects such as the 
appearance of a nontrivial infrared fixed point that goes to infinity at infinite volume or the merg- 
ing of an infrared and an ultraviolet fixed point. We present extensions of the RG flows to the 
complex coupling plane. We discuss the possibility of constructing a continuous beta function 
from the discrete one by using functional conjugation methods. We briefly discuss the relevance 
of these findings for the search of nontrivial fixed points in multiflavor lattice gauge theory mod- 
els. 
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1. Introduction and Background 

With the upcoming experimental effort at the Large Hadron Collider (LHC), there has been a 
lot of activity in the lattice community regarding the existence of non-trivial infrared fixed points 
(IRFP) [1]. Different methods have been used to look for or exclude IRFP and some results are 
controversial. Calculating the discrete j8 function is one of the methods used. In the following, we 
discuss finite volume effects in discrete jS functions for spin models that can be solved by numerical 
methods. 



2. /3 functions and Two Lattice Matching 

In the continuum, the Callan-Symanzik j8 function, /3c5, can be calculated using standai^d 
perturbation theory 

J8cs(«) = -i8ia2-j82a^-i33 + (2.1) 

„2 

where a = f^. The first two coefficients of Eq. 2. 1 are regularization scheme independent. They 
are functions of the number of colors A'^^ and the number of quark flavors Nf [2, 3]: 

l5i = ^{nNc-2Nf) (2.2) 

= 2^ ^^"^^^ " " ^^^^^'^ ^^-^^ 

The theory is asymptotically free if jSi > and a stable IRFP may emerge for /3i > and jSi < 0. 
The IRFP may appear in the large a region where perturbation theory breaks down. Therefore, it 
is essential to locate IRFP nonperturbatively. 

One method to calculate a discrete /3 function nonperturbatively is called the two lattice match- 
ing [5,6]. It is a way of measuring the running of bare couplings based on the fact that all the ob- 
servables will have the same value if the models have the same effective action. One can block-spin 
the system of volume V , where V is the volume in lattice units, n times and calculate an observable 
that will be denoted R{g,V). One can also block-spin the system of smaller volume V /b^ n — I 
times, calculate the observable R{g' ,V /b'^), and tune the new coupling g' until the two observables 
match. If all the observables give the same value for these two coupling configurations g and g', 
then one the two blocked system have the same large distance physics and the same physical cor- 
relation length. In the models discussed below, we only consider one parameter in the coupling 
space. In general, many more couplings are generated during the block-spin transformation. If 
only one relevant parameter exists, all the other irrelevant coupUngs will eventually die out as one 
block-spins many times near the critical point. Therefore, the examples showed below should apply 
to a larger class of system. In Ref. [4], we used the following observable: 

((I;ceBi0x)(LveB2^y)) 

R{l5,r/a^)^' (2.4) 

((I.vEB,0x)(IveB.«^v)))„ 



2 



Volume Effects in Discrete j3 functions 



Yuzhi Liu 



where y is the physical volume of the system, a is the lattice spacing, Bi and B2 denote neighbor 
blocks and D is the dimensionality of the system. The reason we chose this particular rational 
form is that we do not need to calculate the partition function explicitly since it cancels out in 
both denominator and numerator. The matching condition reads R{p,L^) = R{I5', {L/b)'^), where 
b is the scaling factor. The discrete j8 function is defined as AP{P,L^ — )• (L/b)^) = jS — jS'. This 
definition is consistent with the one used in [5, 6]. j(3 denotes either a quantity proportional to g^-^ 
or the inverse temperature and should not be confused with Pes- In order to calculate AjS(j8o), we 
start out from large volume at jSq and tune the coupling at smaller volume {L/b)^ to jSi so that 
R{Po,L^) = R{pi,{L/b)'^). We then match R{Pi,L^) = R{p2, (L/b)'^) and so on. By repeating 
this procedure, we get a sequence of j3's from which we can calculate the discrete j3 function 
defined above. In the following, we will work with Dyson's hierarchical Ising model and the D = 2 
nonlinear 0{N) sigma model in the large Umit. 

3. Dyson's Hierarchical Ising Model 

The Hamiltonian of the Dyson's hierarchical model with 2"'"™ lattice sites is defined as 

^ = -^I(i)"/Wl( I <^^f (3-1) 

«=l xeBi") 

where c controls the interaction strength for different block sizes. We can include the dimension- 
ality D through the relation c/4 = b^^^^ [7], where b = 2^1^ is the scaling factor. The model has 
many nice properties and is an ideal laboratory to test various ideas before one appUes them to the 
more compUcated full QCD case. Here is a list of properties we used [8, 9, 10] 

• For D > 2 and /(«) = 1, the model has a second order phase transition, which is similar to 
the D = 3 regular (nearest neighbor) Ising model. 

• For D <2 and f{n) = 1, the model has no phase transition at finite temperature, which is 
different from the D = 2 regular Ising model. 

• For D = 2 and f{n) = log{n), the model is equivalent to Anderson model and has a Thouless 
effect (discontinuity in the magnetization). 

The cases D= 3, 2, 1.7 with f{n) = 1 have been analyzed thoroughly in [4]. In the following, we 
will consider various dimensions for f{n) = log{n + 1) ( models with /(«) = log{n + 1) have same 
properties as /(«) = log{n) in the infinite volume limit). 

4. How does an IRFP appear/disappear? 

Fixed points correspond to zeros of the j3 function. It has been pointed out that zeros of the jS 
function can disappear in three ways as one or several pai^ameters change [11]. They are schemat- 
ically described in Fig. 1. The first two ways are interchangeable depending on what parameters 
are used (a -h- l/a). For the third way, the fixed points that disappeared can be recovered in the 
complex pai^ameter plane [12]. This third way can be easily seen for the hierarchical model by 
tuning the dimensionality D. Figure 2 shows discrete jS functions for f{n) = log{n + 1), D = 1.9, 
1.994, and 2. The corresponding complex RG flows are shown in Fig. 3. 
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Figure 1: Three mechanisms for the loss of the fixed point(s) described in [11]. 
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Figure 2: Discrete j3 functions for /(«) = log{n + 1) and D=l.9, 1.994, and 2 hierarchical models (dimen- 
sion increases from top to bottom). 



D=1.9,n=4 vs n=5 D=1.994,n=4 vs n=5 D=2,n=4 vs n=5 




Figure 3: Complex RG flows for /(«) = log{n + 1 ) and D=l .9, 1 .994, and 2 hierarchical models. The flows 
are constructed with the two lattice matching method described in the text. 
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D=1.9, f (n) =log (n+l) D=1.994, f (n) =log (n+1 ) D=2, f (n) =log (n+1 ) 
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Figure 4: Discrete j3 functions for /(n) = Zog(« + 1) and D=1.9, 1.994, and 2 hierarchical models. Pseudo 
fixed points may appear for small volumes. Volume dependence is clearly shown in the figure (volume 
increases from bottom to top). 



5. Volume Effects 



In [4], we systematically analyzed complex RG flows for the f{n) = 1 hierarchical model and 
the volume dependence of Fisher's zeros (zeros of the partition function in the complex jS plane). 
Fisher's zeros in lattice gauge models are discussed in [13]. Here we would Uke to illustrate volume 
effects in AjS. For this purpose, we consider the f{n) = log{n + 1) hierarchical model. Figure 
4 shows how discrete j3 functions change with the volume. When the volume is small and the 
dimension parameter D is large enough, pseudo fixed points may appear (Fig. 4 middle and right). 
However as the volume increases, all the zeros of discrete jS functions disappear. This is expected 
since there is no second order phase transition for D = 2 and f{n) = log{n + 1) case. Fixed points 
appearing in smaller volumes are finite volume artifacts. 



6. From Discrete to Continuous Flows 



Up to now, we have discussed properties of discrete /3 functions. If one wants to see how 
the corresponding continuous Callan-Symanzik j8 functions behave, we need to do some interpo- 
lations. In the following, we will construct continuous j3 functions from discrete ones for 0{N) 
models. In [14], complex RG flows and discrete jS functions of 0{N) models have been constructed 
from both two lattice matching and a rescaling method (where the UV cutoff appearing in the di- 
mensionless expression of the bare mass is rescaled). The relation between jS and g is jS g^^ for 
both gauge models and the 0{N) model. By changing the energy scale, j3 — jS' and ^cs — ^ ^'cs' 
we get 

A more general discussion on how to generate continuous flow from step-scaling function has been 
provided in [15]. djS' /dp can be obtained from the discrete j8 function 

1 "^^^ rfi9^ 
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From Eqs. 6. 1 and 6.2, one can easily get Pes from Aj8 by iteration. The only thing one needs 
to fix is the initial starting point for the interaction, i. e. peso- Continuous choice can be obtained 
by requiring that the asymptotic behavior agrees with expansions at small or large coupling. Figure 
5 shows the corresponding continuous j3 functions from both matching and rescaling methods dis- 
cussed in [14]. Colored solid lines are obtained from matching discrete j8 functions with different 
peso- A pseudo fixed point appeared in the small g region and is a finite volume effect. It will go 
away when we further increase the volume. The black dotted line is obtained from the rescaling 
discrete P function. 

Continuous P functions for 2D 0(N) model 
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Figure 5: Continuous j3 functions constructed from both matching and rescaling discrete j3 functions. 



7. Conclusion 

We have successfully constructed discrete P functions for f{n) = log{n + 1) hierarchical Ising 
models for different volumes. It has been clearly demonstrated that volume effects may generate 
pseudo fixed points for systems without phase transition. In order to verify whether the observed 
fixed point(s) is intrinsic or purely finite volume artifact, one needs to do simulations for larger 
volumes. An alternative way is to calculate Fisher's zeros and apply finite size scaling technique. 
The technique of finding Fisher's zeros from the constructed density of state has been developed 
for SU (2) [16] and tested for ?7(1)[17]. The case of SU (3) with various Nf is in progress. 
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